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ABSTRACT 

The lower solar atmosphere consists of partially ionized turbulent plasmas 
harbouring velocity field, magnetic field and current density fluctuations. The 
correlations amongst these small scale fluctuations give rise to large scale 
flows and magnetic fields which decisively affect all transport processes. The 
three fluid system consisting of electrons, ions and neutral particles supports 
nonideal effects such as the Hall effect and the ambipolar diffusion. Here, we 
study magnetic transport by ambipolar diffusion and compare the characteris- 
tic timescales of the laminar and the turbulent ambipolar diffusion processes. 
As expected from a turbulent transport process, the time scale of the turbu- 
lent ambipolar diffusion is found to be smaller by orders of magnitude as com- 
pared with the laminar ambipolar diffusion. The nonlinearity of the laminar 
ambipolar diffusion creates magnetic structures with sharp gradients which 
are amenable to processes such as magnetic reconnection and energy release 
therefrom for heating and flaring of the solar plasma. 



1 INTRODUCTION 

The generation and the transport of magnetic fields, in astrophysical objects in general and in 
the sun in particular, is a subject of great interest and an active area of investigation. Several 
variants of the mean field dynamo have been suggested for the generation and migration of 
the magnetic fiux 

ilsliQ)- In a partially ionized plasma the collisions between the charged 
particles and the neutrals produce additional diffusion mechanisms such as the ambipolar 
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diffusion due to ion-neutral collisions and the magnetic resistivity due to the electron-neutral 
collisions. 

The sun is endowed with a variety of temperature variations resulting from a combi- 
nation of thermal and radiative equilibria and departures therefrom. The thermal and the 
nonthermal nature of processes then translates into a plasma with varying degrees of ion- 
ization. Thus the ionization fraction a = — could vary over several orders of magnitude 

pn 

where pi and pn are respectively the ion and the neutral hydrogen mass densities. Discrete 
structures such as sunspots, prominences and spicules contain plasmas with varying degrees 
of ionization. The support of the neutral component against gravity is a major concern in 
the stability of these structures. Although the ideal magnetohydrodynamics (MHD) is often 
used as a starting point of an investigation, a partially ionized system dominated by the 
charged particle-neutral collisions and the neutral particle dynamics necessitates a 3-fluid 
treatment. The strong charge particle-neutral coupling endows the neutral fluid with some of 
the properties of a conducting fluid. The neutral fluid is thus subjected to the Lorentz force 
along with the usual pressure gradient force. This attribute has been invoked to find the sup- 
port for the neutral component of the partially ionized cold and dense solar prominence (j5|). 
The evolution of the magnetic fields in such a plasma would be affected by the multifiuid 
interactions in general and the ambipolar diffusion in particular!^. The solar magnetic flux, 
generated in the convection zone, has to pass through the partially ionized solar photosphere 
before it can appear high up in the solar corona. This realization is rather recent and is now 
receiving a lot of attention. Arber, Haynes and Leake(0) has emphasized the profound effects 
on the temperature and the current structure of the overlying chromosphere and the corona 
that the inclusion of the neutral medium can produce. Vogler and Schiisslerlsl) have invoked 
dynamo mechanism in a region a few hundred kilometers below and above the visible solar 
surface to account for the stronger horizontal magnetic fields. The issue of possible discon- 
nection between the sub-surface and the surface solar magnetic field, recently emphasized 
by Schiisslerjoj), may have some bearing on the neglect of the neutral fiuid-plasma coupling 
in the flux transport on the solar photosphere. It is clear that the transport and or genera- 
tion processes of the magnetic field on the solar photosphere must be studied in a partially 
ionized plasma. The magnetic transport would occur through the large scale flows as well 
as the turbulent fluctuations of the velocity field, the magnetic field and the current density 
with which the photosphere is well endowed. Recently Krishan and Gangadhara(llOl) have 
initiated the study of mean-field dynamo in a partially ionized plasma. 
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For a low degree of ionization one can define a weakly ionized plasma by the condition (j 111) 
that the electron-neutral collision frequency z/e„ ~ 10~-^^nn.( ^^^"^ )°'^ is much larger than the 
electron- ion collision frequency z/gj ~ 6 x 10~'^'^niAZ'^{KBT)~^'^ . This translates into the 



ionization fraction Up/un < 5 x 10^ T where n's are the particle densities and T is the 
temperature in Kelvin. A major part of the solar photosphere (j5|) qualifies as a weakly ionized 
plasma (1 13l ). 

In this paper, we investigate the magnetic transport by the flow generated by the mag- 
netic field itself through the ambipolar diffusion in a self consistent manner. We develop a 
three fluid framework in section two. The turbulent electromotive force and the large scale 
field dynamics is established in section 3. The advection of the large scale magnetic field, in 
the kinematic limit Il4 ) is studied and the timescales of advection by the laminar and the 
turbulent ambipolar diffusion are presented in section 4. We end the paper with a section 
on conclusion. 



2 THREE-COMPONENT MAGNETOFLUID 



We begin with the three component weakly ionized plasma consisting of electrons (e), ions 
(i) of uniform mass density pi and neutral particles (n) of uniform mass density p„. The 
equation of motion of the electrons can be written as: 



dt 



E + 



V,xB 



-Vpe- 



(8.1) 



where the electron-ion collisions have been neglected since the ionized component is of low 
density. On neglecting the electron inertial force, the electric field E is found to be: 
V,x B Vpc 



E 



V^) . (8.2) 

c eue e 

This gives us Ohm's law. For time scales of interest larger than the ion-neutral collision time 

scale the ion dynamics can be ignored. The ion force balance then becomes: 

Vi X B 



= —Vpi + eui 



E 



^inPiiy, - Vn) , (8.3) 

where v\ri is the ion-neutral collision frequency, and the ion-electron collisions have been 
neglected for the low density ionized component. Substituting for E from Eq. (2) we find 



the relative velocity between the ions and the neutrals: 



lA) 
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where 

J = ene(Vi - Ve) . 

The equation of motion of the neutral fluid is: 



(8.5) 



Pn 



dt 



-Vpn - «^niPn(l^n " ^i) - 
'^nePn(T'^n ^e) ; 



(8.6) 



where the viscosity of the neutral fluid has been neglected. Substituting for — M from 
Eq. (4), and using z/jnPi = Vnipn we find: 



Pn 



dV 



dt 



- + (Vn- V)V„ 



J X B 



(8.7) 



where p = Pn + Pi + Pe- Observe that the neutral fluid is subjected to the Lorentz force as a 
result of the strong ion-neutral coupling due to their collisions. 
Consider Faraday's law of induction: 



dB 
~dt 



-cV X E 



By substituting for the electric field from Eq. (2), we get 



dB 
~dt 



V X (Ve X S) + 7]V'^B 



18) 



19) 



where the pressure gradient terms have been dropped for the incompressible case with con- 
stant temperature. Here r] = meZ/enC^/(4vre^ne) is the electrical resistivity predominantly due 
to electron-neutral collisions. Using the construction 



X B = [V„ - (V„ - Vi) - (Vi - Ve)] X B 



i.lO) 



and substituting for the relative velocity of the ion and the neutral fluid from Eq. (4), Eq. (9) 
becomes: 



dB 



V X 



J JxB\ „ 

Vn + X B 

eUe CZ/inPi / 



+ T]V^B 



(8.11) 



One can easily identify the Hall term {J/eUe), and the ambipolar diffusion term (J x -B)(j7|) 
. The Hall term is much larger than the ambipolar term for large neutral particle densities 
or for z/in ^ ujci where uJd is the ion cyclotron frequency. In this system the magnetic field is 
not frozen to any of the fluids. Equations (7) and (11) along with the mass conservation 



V ■ = 

form the basis of our investigation. 



(8.12) 
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3 TURBULENCE IN THREE-COMPONENT MAGNETOFLUID 

One way of studying turbulence in a system is by splitting the physical quantities into small 
scale and large scale parts. The mean-field dynamo(.4) is a consequence of such a study. We 
will tread the same path. The magnetic induction equation (11) is written as: 

V X X + r/V^S , (8.13) 



dt 
where 

Ve = V„ + Vh + Vah. (8.14) 
with 

Vh = -— (8.15) 
as the Hall velocity and 

VAm= 8.16 

could be called the ambipolar velocity. Following the standard procedure (jj) the velocity Ve 
and the magnetic field B are split into their average large scale parts and the fiuctuating 
small scale parts as: 

Ve = V^+n, (8.17) 
B = B + B' (8.18) 

such that 

V^ = 0, W=0. (8.19) 

In the kinematic dynamo the magnetic induction equation is solved for large and small scale 
fields. Substituting Eqs. (17) and (18) into the induction equation (11), we find, in the first 
order smoothing approximation, 

n^K-^ + ^H-^ (8.20) 

eUe Cl/inPi CZ/inPi 

and the mean fiow is found to be: 



J J X B J' X B' , , 

Ve = + + . 8.21 

erie CZ/inPi CZ/inPi 

The turbulent electromotive force S is a function of the mean magnetic induction B and 
mean quantities formed from the fiuctuations, and is expressed as: 



E = V'eX B' = aB - X B , (8.22) 
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where 



a — 



Here 



-n • (V X n) 



3 

CKv + CKh + CKAm ■ (8.23) 



«v = -^K • Q', (8.24) 

is the measure of the average kinetic hehcity of the neutral fluid in the turbulence possessing 

correlations over time Tcor and 
2r„ 



an = =^J' ■ n', (8.25) 

represents the contribution of the Hall effect. The coupling of the charged components with 
the neutral fluid is clearly manifest through the possible correlation between the current 
density fluctuations and the vorticity fluctuations of the neutral fluid 17^ = V x V^. The 
ambipolar term gives rise to 

"Am = OiA B , (8.26) 

with 

2re, 



aA = ^-^J' X , (8.27) 

3cpii/in 

as the contribution from the ambipolar diffusion with its essential nonlinear character man- 
ifest through its dependence on the average magnetic induction. One also observes that the 
Hall alpha (Eq. 25) requires a component of the fluctuating current density along the fluc- 
tuating vorticity of the neutral fluid whereas the ambipolar effect (Eq. 27) thrives on the 
component of the fluctuating current density perpendicular to the fluctuating vorticity. The 
turbulent dissipation is given by 



P = -yV'i ^f3, + Pu + /3Am (8.28) 
with 



/3v = (8.29) 

as the measure of the average turbulent kinetic energy of the neutral fluid in the turbulence 

possessing correlations over time Tcor and 
2r, 



Ph = -TT^ J' • K (8.30) 

SeUe 

represents the contribution of the Hall effect. The coupling of the charged components with 
the neutral fluid is clearly manifest through the possible correlation between the current 
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density fluctuations and the velocity fluctuations of the neutral fluid. The ambipolar term 
furnishes 

PAm^(3A-B, (8.31) 
/3a = J^TV^ (8.32) 

with its essential nonlinear character manifest through its dependence on the average mag- 
netic induction. One also observes that the Hall /9h requires a component of the current 
density fluctuations along the velocity fluctuations of the neutral fluid whereas the am- 
bipolar effect thrives on the component of the current density fluctuations perpendicular to 
the velocity fluctuations. We have used rigid or perfectly conducting boundary conditions 
(all surface contributions vanish) while determining the averages. The magnetic transport 
equation becomes: 
dB 



dt 



V X 



VeX B + aB - (3V X B + rjV'B . (8.33) 



4 MAGNETIC TRANSPORT BY TURBULENT AMBIPOLAR 
DIFFUSION 

In order to highlight the effect of the ambipolar diffusion on the magnetic transport, We 
ignore all other effects so that the large scale neutral fluid flow Vn = 0, the Hall flow 
Vh = — ^ = 0,q; = 0,/3 = and 77 = 0. The transport equation reduces to 

dB 



dt 



= V X 



Jx B — J' xB' — 
X B-\ X B 



(8.34) 



The flrst term on the right hand side represents the laminar ambipolar diffusion. The second 
term represents the turbulent ambipolar diffusion or transport of the magnetic induction. 
We estimate the second term as follows: 



B'.VB' 



J' xB'^ + c — . (8.35) 

87r 47r 

The left hand side is the average Lorentz force due to magnetic fluctuations and we assume 
that it is nonzero. On the right side, the first term is symptomatic of the turbulent mag- 
netic pressure and the second is the contribution to the turbulent curvature. If magnetic 
fluctuations are introduced in a fluid initially at rest, the Lorentz force will drive motion 
in the fluid. In an incompressible fluid the direction of the flow is along the predominant 
component of the fluctuating fleld. This arises through the curvature term B'.VB' since 
only the solenoidal part of the Lorentz force can drive flows in an incompressible fluid or in 
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a compressible fluid under the Boussinesq approximation as discussed by Ogilvie(Q). The 
gradient part of the Lorentz force is compensated by the pressure gradient. Here we shall 
make an approximate estimate of this term as: 



c 



J' xB' ^ —^B'^ . (8.36) 



where Xcor is the correlation length of the magnetic fluctuations. We follow ('4') in order to 
estimate the mean square of the magnetic fluctuations B''^. In the high conductivity limit 
(r] ^ 0), the fluctuation B' is produced by the interaction of the turbulent velocity field V'^ 
with the mean magnetic field B over a time scale Tcor where Tcot is the correlation timescale 
of the fluctuating fields. Thus the elemental field B'^i is given by 

B'^, - . (8.37) 

After a time interval Tcor the turbulent field V'^ changes nearly completely. This new realiza- 
tion of the turbulent V'^ interacts with B to produce another elemental field B'^^ uncorrelated 
with the previous field element. Each field element survives over the dissipation time scale 
Tdis ~ much larger than the correlation time Tcor- Thus the total fluctuating field B' is 
the incoherent sum of these elemental fields B'^i. The number n of these elemental fields is 
given by 

n ^ ^ . (8.38) 
Thus, to the order of magnitude, B' is found to be 

B' » VSBrf 



^ '°'B . (8.39) 



V 



And 

fi^ ^ RmSB'^ . (8.40) 
where 

= Yhh^^ s = . (8.41) 

are respectively the magnetic Reynolds number and the Strouhal number. We can now 
estimate the timescales of the laminar ambipolar diffusion, Tla and the turbulent ambipolar 
diffusion, Tta from the magnetic transport equation (Eq.34) as: 

AnL'^PiUin AnLXcorPi^in 

Tla =2 , Tta = • (8.42) 

where L is the characteristic spatial scale associated with the mean magnetic field B. Sub- 
stituting from Eq.(40), we find the ratio of the two timescales to be: 
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^ = (RMSr^ . (8.43) 

Now the magnetic Reynolds number Rm >> 1, the Strouhal number S < 1 and the corre- 
lation length Xcor << L] therefore the turbulent ambipolar diffusion timescale Tta can be 
much smaller than the laminar ambipolar diffusion timescale Tl^. This is exactly what is 
expected of a turbulent transport process. The reason for the greater efficiency of a turbulent 
transport over a laminar transport is that in a turbulent process a physical attribute on a 
large scale is transferred to a small scale and it is much easier for a small scale quantity 
to diffuse than for a large scale one. For the typical values of the physical parameters at a 
height of 500Km on the solar atmosphere we find T^a ~ lO^sec and Tta could be smaller 
by several orders of magnitude. 



5 CONCLUSION 

It is found that the ambipolar diffusion is an important effect in the partially ionized part 
of the solar atmosphere. The laminar ambipolar transport is known to create magnetic 
structures ({gI) with steep gradients(? ). The turbulent part of the ambipolar diffusion could 
transport these structures on a timescale which could be shorter by orders of magnitude 
than all other processes. Thus the formation and migration of magnetic structures to higher 
up in the solar atmosphere could be the result of ambipolar diffusion. 
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